COVID-19 epidemic modelling for policy decision support in Victoria, Australia 2020-2021.

BACKGROUND: Policy responses to COVID-19 in Victoria, Australia over 2020-2021 have been supported by evidence generated through mathematical modelling. This study describes the design, key findings, and process for policy translation of a series of modelling studies conducted for the Victorian Department of Health COVID-19 response team during this period. METHODS: An agent-based model, Covasim, was used to simulate the impact of policy interventions on COVID-19 outbreaks and epidemic waves. The model was continually adapted to enable scenario analysis of settings or policies being considered at the time (e.g. elimination of community transmission versus disease control). Model scenarios were co-designed with government, to fill evidence gaps prior to key decisions. RESULTS: Understanding outbreak risk following incursions was critical to eliminating community COVID-19 transmission. Analyses showed risk depended on whether the first detected case was the index case, a primary contact of the index case, or a 'mystery case'. There were benefits of early lockdown on first case detection and gradual easing of restrictions to minimise resurgence risk from undetected cases. As vaccination coverage increased and the focus shifted to controlling rather than eliminating community transmission, understanding health system demand was critical. Analyses showed that vaccines alone could not protect health systems and need to be complemented with other public health measures. CONCLUSIONS: Model evidence offered the greatest value when decisions needed to be made pre-emptively, or for questions that could not be answered with empiric data and data analysis alone. Co-designing scenarios with policy-makers ensured relevance and increased policy translation.

Refining Reproduction Number Estimates to Account for Unobserved Generations of Infection in Emerging Epidemics.

BACKGROUND: Estimating the transmissibility of infectious diseases is key to inform situational awareness and for response planning. Several methods tend to overestimate the basic (R0) and effective (Rt) reproduction numbers during the initial phases of an epidemic. In this work we explore the impact of incomplete observations and underreporting of the first generations of infections during the initial epidemic phase. METHODS: We propose a debiasing procedure that utilizes a linear exponential growth model to infer unobserved initial generations of infections and apply it to EpiEstim. We assess the performance of our adjustment using simulated data, considering different levels of transmissibility and reporting rates. We also apply the proposed correction to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) incidence data reported in Italy, Sweden, the United Kingdom, and the United States. RESULTS: In all simulation scenarios, our adjustment outperforms the original EpiEstim method. The proposed correction reduces the systematic bias, and the quantification of uncertainty is more precise, as better coverage of the true R0 values is achieved with tighter credible intervals. When applied to real-world data, the proposed adjustment produces basic reproduction number estimates that closely match the estimates obtained in other studies while making use of a minimal amount of data. CONCLUSIONS: The proposed adjustment refines the reproduction number estimates obtained with the current EpiEstim implementation by producing improved, more precise estimates earlier than with the original method. This has relevant public health implications.

A Comparative Analysis of Statistical Methods to Estimate the Reproduction Number in Emerging Epidemics, With Implications for the Current Coronavirus Disease 2019 (COVID-19) Pandemic.

BACKGROUND: As the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic continues its rapid global spread, quantification of local transmission patterns has been, and will continue to be, critical for guiding the pandemic response. Understanding the accuracy and limitations of statistical methods to estimate the basic reproduction number, R0, in the context of emerging epidemics is therefore vital to ensure appropriate interpretation of results and the subsequent implications for control efforts. METHODS: Using simulated epidemic data, we assess the performance of 7 commonly used statistical methods to estimate R0 as they would be applied in a real-time outbreak analysis scenario: fitting to an increasing number of data points over time and with varying levels of random noise in the data. Method comparison was also conducted on empirical outbreak data, using Zika surveillance data from the 2015-2016 epidemic in Latin America and the Caribbean. RESULTS: We find that most methods considered here frequently overestimate R0 in the early stages of epidemic growth on simulated data, the magnitude of which decreases when fitted to an increasing number of time points. This trend of decreasing bias over time can easily lead to incorrect conclusions about the course of the epidemic or the need for control efforts. CONCLUSIONS: We show that true changes in pathogen transmissibility can be difficult to disentangle from changes in methodological accuracy and precision in the early stages of epidemic growth, particularly for data with significant over-dispersion. As localized epidemics of SARS-CoV-2 take hold around the globe, awareness of this trend will be important for appropriately cautious interpretation of results and subsequent guidance for control efforts.

Reconstructing the early global dynamics of under-ascertained COVID-19 cases and infections.

BACKGROUND: Asymptomatic or subclinical SARS-CoV-2 infections are often unreported, which means that confirmed case counts may not accurately reflect underlying epidemic dynamics. Understanding the level of ascertainment (the ratio of confirmed symptomatic cases to the true number of symptomatic individuals) and undetected epidemic progression is crucial to informing COVID-19 response planning, including the introduction and relaxation of control measures. Estimating case ascertainment over time allows for accurate estimates of specific outcomes such as seroprevalence, which is essential for planning control measures. METHODS: Using reported data on COVID-19 cases and fatalities globally, we estimated the proportion of symptomatic cases (i.e. any person with any of fever ≥ 37.5 °C, cough, shortness of breath, sudden onset of anosmia, ageusia or dysgeusia illness) that were reported in 210 countries and territories, given those countries had experienced more than ten deaths. We used published estimates of the baseline case fatality ratio (CFR), which was adjusted for delays and under-ascertainment, then calculated the ratio of this baseline CFR to an estimated local delay-adjusted CFR to estimate the level of under-ascertainment in a particular location. We then fit a Bayesian Gaussian process model to estimate the temporal pattern of under-ascertainment. RESULTS: Based on reported cases and deaths, we estimated that, during March 2020, the median percentage of symptomatic cases detected across the 84 countries which experienced more than ten deaths ranged from 2.4% (Bangladesh) to 100% (Chile). Across the ten countries with the highest number of total confirmed cases as of 6 July 2020, we estimated that the peak number of symptomatic cases ranged from 1.4 times (Chile) to 18 times (France) larger than reported. Comparing our model with national and regional seroprevalence data where available, we find that our estimates are consistent with observed values. Finally, we estimated seroprevalence for each country. As of 7 June, our seroprevalence estimates range from 0% (many countries) to 13% (95% CrI 5.6-24%) (Belgium). CONCLUSIONS: We found substantial under-ascertainment of symptomatic cases, particularly at the peak of the first wave of the SARS-CoV-2 pandemic, in many countries. Reported case counts will therefore likely underestimate the rate of outbreak growth initially and underestimate the decline in the later stages of an epidemic. Although there was considerable under-reporting in many locations, our estimates were consistent with emerging serological data, suggesting that the proportion of each country's population infected with SARS-CoV-2 worldwide is generally low.

SBDiEM: A new mathematical model of infectious disease dynamics.

A worldwide multi-scale interplay among a plethora of factors, ranging from micro-pathogens and individual or population interactions to macro-scale environmental, socio-economic and demographic conditions, entails the development of highly sophisticated mathematical models for robust representation of the contagious disease dynamics that would lead to the improvement of current outbreak control strategies and vaccination and prevention policies. Due to the complexity of the underlying interactions, both deterministic and stochastic epidemiological models are built upon incomplete information regarding the infectious network. Hence, rigorous mathematical epidemiology models can be utilized to combat epidemic outbreaks. We introduce a new spatiotemporal approach (SBDiEM) for modeling, forecasting and nowcasting infectious dynamics, particularly in light of recent efforts to establish a global surveillance network for combating pandemics with the use of artificial intelligence. This model can be adjusted to describe past outbreaks as well as COVID-19. Our novel methodology may have important implications for national health systems, international stakeholders and policy makers.